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Abstract 

We study the instabilities of quark matter in the framework of a gener- 
alized Nambu-Jona-Lasinio model, in order to explore possible competition 
between three-quark clustering to form nucleons and diquark formation lead- 
ing to color superconductivity. Nucleon and A solutions are obtained for the 
relativistic Faddeev equation at finite density and their binding energies are 
compared with those for the scalar and axial- vector diquarks found from the 
Bethe-Salpeter equation. In a model with interactions in both scalar and 
axial diquark channels, bound nucleons exist up to nuclear matter density. 
However, except at densities below about a quarter of that of nuclear matter, 
we find that scalar diquark formation is energetically favored. This raises 
the question of whether a realistic phase diagram of baryonic matter can be 
obtained from any model which does not incorporate color confinement. 
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I. INTRODUCTION 



In contrast to QCD at finite temperature, rather little is known about QCD at finite 
density. Technical difficulties (such as the complex fermionic determinant at finite chemical 
potential) make lattice Monte-Carlo simulations very difficult. Even though some techniques 
are being developed to overcome these problems (for example, the Glasgow method M or the 
technique of imaginary chemical potential 0), they are not yet able to provide unambiguous 
results. 

However, models of QCD seem to indicate a rich phase structure in high-density quark 
matter. In particular, much attention has recently been devoted to so-called "color su- 
perconductivity" [HIJ: an arbitrarily weak attractive force makes the Fermi sea of quarks 
unstable at high density with respect to diquark formation and induces Cooper pairing (di- 
quark condensation). Although this phenomenon had been studied earlier ||||[8|], the large 
magnitude of the superconducting gap found in the more recent studies 0,f|] (more than 
100 MeV) suggested that this was much more important than had been thought previously 
and has generated an extensive literature. In Refs. |||| an instanton model was used to 
calculate the gap at finite density and zero temperature. Berges and Rajagopal extended 
this work and calculated the phase diagram of strongly interacting matter as a function of 
temperature and baryon-number density in the same model. Identical results were found 
in the Nambu-Jona-Lasinio model ||. The existence of such a color-superconducting gap 
could have important consequences for the physics of neutron stars or even for heavy ion 
collisions [10|,11|. (For a review of the field, see: [p~2]l .) 

Previous studies of color superconductivity have focused on instabilities of the quark 
Fermi sea with respect to diquarks only. Of course, it is known that at lower densities 
(of the order of nuclear matter density) three-quark clusters — nucleons — are the dominant 
degrees of freedom. Here we address the question of the possibility of a competition between 
diquark condensation and three-quark clustering at finite density. To answer such a question 
fully would require a three-particle generalization of the BCS treatment, a complicated task. 
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As a first step towards this goal, we look for instabilities of the quark Fermi sea with respect 
to three-quark clustering by studying the nucleon binding energy in quark matter with finite 
density. A bound nucleon at finite density would be a signal of instability with respect to 
clustering (in the same way that a bound diquark at finite density is a signal of instability 
with respect to diquark condensation). A comparison of the magnitudes of the binding 
energies of the diquark and the nucleon can give some idea of the relevance of these degrees 
of freedom. Very recently, Beyer et al. have studied a similar clustering problem for three 
nucleons in nuclear matter fl3| . 



To perform this study we solve the Bethe-Salpeter equation for the diquarks and the Fad- 
deev equation for the nucleon (and also, for completeness, for the A) within the framework 
of a Nambu-Jona-Lasinio (NJL) model. The use of a separable interaction simplifies the 
treatment of the Faddeev equation considerably. The form of the equation reduces to that 
of a Bethe-Salpeter equation describing the interaction between a quark and a diquark. At 
zero density, several groups have used this Faddeev approach to study baryons in the NJL 
model (for example, Refs. fln)-|iT|). More recent papers ^U] have extended this treatment to 



incorporate a mechanism for confinement. In this work we will mainly follow the formalism 
developed by Ishii, Bentz and Yazaki generalizing it to finite density. In this type of 



approach, it is crucial to include the axial-vector diquark channel in addition to the scalar 
channel, as otherwise the nucleon is very weakly bound. In previous preliminary studies, 
keeping only the scalar channel, we have found that binding of a nucleon in matter is only 
possible at very low densities, less than 10% of nuclear matter density ||21|| . 



The paper is organized as follows: in the next section, we briefly describe the NJL model 
and its application to quark-quark interaction. The parameters used in our study are also 
discussed. In Sec. 3 the Bethe-Salpeter equations in the scalar and axial-vector diquark 
channels are solved. The form of the Faddeev equation for the nucleon is presented in 
Sec. 4. In Sec. 5, we describe the numerical techniques used and in Sec. 6 we present our 
results for the nucleon and the A at finite density. Finally, we draw some conclusions in 
Sec. 7. 



II. THE MODEL 



The NJL model provides a simple implementation of dynamically broken chiral symme- 
try, based on a two-body contact interaction f2"2|| . In spite of the fact that the model does 
not incorporate confinement, it has been successfully applied to the description of mesonic 
properties at low energy. (For reviews, see Refs. p3jp4| .) 

The model Lagrangian has the form 

C = - m)V + Ci, (1) 

where Cj is the interaction Lagrangian. In the present work we consider only the chiral 
limit, setting the current quark mass m to zero. Several versions of the NJL interaction can 



be found in the literature. The original version [£2] is 

d = g[(W) 2 + @i l5 T^) 2 ]. (2) 
One can also work with a color-current interaction, 

£/ = "SE(?74 A ^) 2 , (3) 

c 

where A c (c = 1, 8) are the usual Gell-Mann matrices. 

Whatever version of the model is chosen, a Fierz transformation should be performed 
in order to antisymmetrize the interaction Lagrangian. This allows Ci to be brought into 
a form where the interaction strength in a particular channel can be read off directly from 
its coefficient in the Lagrangian. For the qq channel, one just rewrites Cj into the form 
£-i,qq = \{Ci + £i,f) where Cj } p is the Fierz transformed form of £/. We shall need here 
only the scalar and pseudoscalar terms of the qq interaction. These have the same form as 
Eq. (0), but with a coupling constant g n which is related to the original coupling constant g 
of the Lagrangian by a coefficient given by the Fierz transformation. For example, one has 
g n /g = 13/12 for the model defined by Eq. (g) and g v /g = -2/9 for Eq. (|). 

To study the nucleon we also need to rewrite the interaction Lagrangian in the form of 
a qq interaction. This is done by a Fierz transformation to the qq channels, which allows 



the interaction to be expressed as a sum of terms of the form (ipAip )(ip T Bip), where the 
matrices A and B are overall antisymmetric in Dirac, isospin and color indices. (We use the 
Dirac representation for the 7-matrices and follow the conventions of Itzykson and Zuber 
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As our three-quark state must be a color singlet, the diquark channels of interest are 
color anti-triplet. For a local interaction the relevant channels are the scalar (0 + , T = 0) and 
axial- vector (1 + , T = 1) ones. These are also the channels in which more realistic interactions 
(including, for example, one-gluon exchange) are expected to be most attractive. Explicitly, 
we have 

C s = g s £ (^7 5 Cr 2 /?> T ) (^C" , (4) 

a 

for the scalar channel and 

£a = 0aE (HWf) (> T C- Vr 2 r 4 /3>) , (5) 

i,a 

for the axial- vector one. The matrices [3 a = y| A a for a = 2, 5, 7 project onto the color 3 
channel and C = 27270 is the charge conjugation matrix. The coupling strengths g s and g A 
are again related to the original g by a coefficient given by the Fierz transformation. In the 
following we do not choose a specific version of the NJL Lagrangian but instead treat the 
physical couplings g n ,gs,9A as independent parameters. 

The gap equation for the constituent quark mass M reads 

M = 2ig v f^ i Tr[S(k)] 1 (6) 



where the quark propagator is 



The integral (||) diverges and so has to be regularized. There are various regularization 
schemes at our disposal: Pauli-Villars, proper-time, and 3- or 4-momentum cut-off. In this 
work we use a sharp cut-off A on the 3-momentum to regularize the loop integrals, since 



TABLE I. Values of the parameters g n and A for the two values we use for the constituent 
quark mass M. For each M, the qq coupling ratios gs/g-K and gA/dn are determined by fitting the 
N and A masses. In each case we also list, in parentheses, the minimal value of the coupling ratio 
required to produce bound diquarks. 



M 


9- 


A 


9s 1 9* (min.) 


9a/9tt (min.) 


(GeV) 


(GeV~ 2 ) 


(GeV) 






0.45 


7.914 


0.579 


0.72 (0.33) 


0.46 (0.37) 


0.5 


8.634 


0.573 


0.72 (0.28) 


0.53 (0.31) 



this is conveniently applied to systems of finite density. Although a 3-momentum cut-off 
is not Lorentz invariant, this is less relevant at finite density where there is a much larger, 
physical, breaking of Lorentz invariance due to the presence of a quark Fermi sea. In any 
case, we shall show that physical observables depend only weakly on the choice of regulator. 

The two parameters g n and A can be fitted to a given value of the constituent quark 
mass and to the pion decay constant f n = 93 MeV. This last quantity is evaluated from 

In the following calculations we use two different values for the constituent mass M = 450 
MeV and M = 500 MeV, both of which have been chosen to be large enough that the mass 
of the A lies below the three-quark threshold. In Table | we list the values of g n and A for 
the corresponding constituent masses. As one can see, the values for the cut-off are relatively 
low and decrease with increasing constituent mass. One cannot push the constituent mass 
to higher values as otherwise the cut-off would approach (or, worse, become smaller than) 
the quark mass. 

The effects of finite density on the constituent mass are taken into account by introducing 
the Fermi momentum kp as a lower cut-off on the integral in the gap equation (|j). In Figure 
|l] we show the evolution of the constituent mass as a function of the Fermi momentum for 
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0.5 




k F (GeV) 

FIG. 1. Evolution of the constituent quark mass M with the Fermi momentum kp for 
M(k F = 0) = 450 MeV (continuous curve) and M(k F = 0) = 500 MeV (dashed). 

the two sets of parameters corresponding to M = 450 and 500 MeV. For these parameters, 
the restoration of chiral symmetry occurs at kp = 359 MeV and kp = 372 MeV respectively 
(For comparison, nuclear matter density corresponds to a quark Fermi momentum of kp = 
[§tt 2 pb] 1/3 = 270 MeV). 

III. THE TWO-BODY T-MATRIX 

The diquark T-matrix is an essential building block for the Faddeev equation. It is 
obtained by solving the Bethe-Salpeter equation in the ladder approximation, 

T a(3 ^{k) = K a ^ 5 + - V / -^-K a ^ Xe S xx '(k + q)S et '(-q)T x ' t ''^{k). (9) 

2 \,e,\>,e> J ( 27F ) 

In the scalar diquark channel, the interaction kernel K is 

Kf« s = Aig s ^sCr^T^C-^r^y 6 . (10) 

a 

This interaction is momentum-independent and so Eq. (^) can be solved easily to get the 
diquark T-matrix in the scalar channel jl8fl : 
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(11) 



with 



f s (k) 



l + 2g s U s (k 2 Y 



(12) 



and 



n s (A; 2 ) =Gif ^-Jr D [ l5 S(qh 5 S(k + q)}. 



(27T) 



(13) 



If the T-matrix fll2|) has a pole, this gives us the mass of the bound scalar diquark. Note 
that, if one replaces the coupling g s by g n , the denominator of (|12"D is the same as that 
in the pion qq channel |^8| . This means that for g s = g n the scalar diquark and pion are 
degenerate, and have zero mass in the chiral limit. Since diquarks do not condense in the 
vacuum, this puts an upper limit to the choice of the scalar coupling g s . 
In the axial-vector channel, the kernel is 



(14) 



and the solution of (ft) can be shown to be 



(15) 



with 



7f (fc) = 



1 + 2(? A n A , T (P) 1 + 2g A Tl A , L (k 2 ) 
Here, the axial polarization tensor, 

d A q 



has been decomposed in the form 



n^(k) = n 4 , T (fc 2 ) _ ** + n 4 , L (^ 



k 2 



k 2 



(16) 



(17) 



(18) 
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Again, a bound axial- vector diquark corresponds to a pole in the T-matrix (0). The 
minimum values of the coupling ratios g 3 /g^ and g A J g n required to get bound diquarks at 
zero density can be found in Table I. 

Note that when the loop integrals are regulated with a simple cut-off on either the 3- 
or 4-momentum, the longitudinal polarizability in this channel, II 4 L (A; 2 ) does not vanish, 
in contrast to the hybrid method involving dimensional regularization used by Ishii et al. 
18| , |26|| . However, since there is no conserved current coupled to these states, this does not 



violate any physical symmetries. 

In the presence of quark matter, manifest Lorentz invariance is broken and the structure 
of the polarization tensor Xl^ik) becomes more complicated. For a diquark momentum in 
the z-direction, k = (k°, 0, 0, k 3 ), it can be written 

/ 



u3\ 



.4 











B 











B 


Dk°k 3 









Dk°k 3 



B + C(k 3 ) 2 



(19) 



The use of a 3-momentum cut-off does lead to deviations from the Lorentz covariant structure 
shown in (1181) even in the vacuum case, but we have checked that these are small. 



The results for diquarks at zero density are qualitatively similar to those in Ref. [18 
despite a different choice of regulator and the use of a nonzero current-quark mass in that 
work. For a constituent quark mass of M = 400 MeV, as used in that work, we find that 
the minimum value of g s / 'g n for diquark binding is 0.4 compared with 0.33 in Ref. |l8j. For 
gsl g-n = 0.6 and g s j ' g n = 0.8 we get scalar diquark masses of 699 and 507 MeV respectively, 



compared with 627 and 446 MeV (cf. Table 2 of Ref. ]T8|). To bind the axial- vector diquark, 
we find a minimum coupling strength g A /g-K = 0.46. The corresponding value given in 
Ref. is 2.0, but this should in fact be divided by a factor of 4 p6| , and so we again have 
qualitative agreement with that work. 
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p+P/2 



p'+P/2 



p"+P/2 



= -p-p + -p-p 

-p+P/2 -p'+P/2 — .p"+P/2 

FIG. 2. Diagrammatic representation of the Faddeev equation (|20|). 
IV. THE FADDEEV EQUATION 

Because of the separability of the NJL interaction, the ladder approximation to the 
Faddeev equation for the three-body system can be reduced to an effective two-body Bethe- 
Salpeter equation. This can be thought of as describing the interaction between a quark 
and a diquark, although it is not necessary that the diquark be bound. In our derivation 
of this equation, we have followed the procedure described in We reproduce here only 



the main steps of this derivation; the details can be found in the original papers fllq , |18l . In 
the case of a purely scalar qq interaction, a thorough discussion of the Faddeev equation can 



also be found in Ref. [fT^] . 

We denote the scattering amplitude of a quark on a diquark by '. The indices a, b label 
the diquark and, following the convention of Ref. [DJ they take the values 5 (for the scalar 
diquark) and 0, 3, +1 and —1 (for the components of the axial-vector diquark in a spherical 
basis). These diquark indices will be written as subscripts or superscripts to indicate that 
the components of the axial diquark are covariant or contravariant respectively. The Dirac 
indices a,/3 label the quark, taking the values 1 to 4. The amplitude obeys the integral 
equation, 

x£tf,p) = Z$(j/,p) + £ J ^^(p',/)^(P/2 + /)f cd (P/2-/)Xf(/,p), 



y,c,5,d ' 



(20) 



corresponding to the diagram shown in Fig. |2|. The piece of the kernel containing the 
propagator of the exchanged quark is 



(21) 



7,<5 
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where Q and Q are the two-body vertex functions, already used in Eqs. fllOD and ([14]). The 
propagator of the spectator quark is S(P/2+p"), and f cd (P/2—p") is the two-body amplitude 
for the two interacting quarks. The separable nature of the two-body interaction means that 
this amplitude can be thought of as a diquark propagator, but one should remember that it 
describes the propagation of all two-quark states, not just the bound states (if they exist). 
It can be written as 



-cd 



'r. o \ 







(22) 



where f s and f^' are the scalar and axial-vector diquark "propagators" given by Eqs. flT2"D 
and QI5D respectively. 

To find bound states of the three-body system, we solve the homogeneous version of 



equation ( P0|) for the effective two-body vertex function of a quark and a diquark. This 
vertex function, X", is related to the amplitude by 

x:(p')x^ P ) 



as 



P l 



M 



Bi 



P 2 - M 2 B + ie 

where Mb is the mass of the bound baryon. From ( |20| ) one obtains the equation 
X a(p)= E j ^Z:2{p,pi)S^\P/2+p'^ d {P/2-pi)X 5 d {pi), 



(23) 



(24) 



wJ / (2tt) 
for the vertex function. 

This equation needs to be projected onto states of definite color, spin and isospin. Pro- 
jecting the kernel Z" c 7 onto a color-singlet state gives 



z:?(p, P ') = -3 

for the isospin- 1 channel, and 



1 l^S (p + p')75 V&YnS (p + j/)j s ^ 



V3~f 5 S (p + p')7m' -7a* 5 (P + P'hh 



(25) 



/ 



07 



z:2(p,p>) = -6( 1 ,s(p+p') 1 ,,)^ 



(26) 



for the isospin- § channel. 
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Next, a projection onto good spin and parity must be carried out. Like Ishii et al. [18 
we use the helicity formalism of Jacob and Wick [Q , constructing first a basis of states with 
definite helicity by acting with a rotation operator on helicity eigenstates whose momentum 
p lies along the z-axis: 



| P H,«,a) = K(u)\f>,a,a) = £ S a ' a (uj)R a ' a (u)\p(uj), a', a), 



(27) 



where p(u) = R(u)p lies in a general direction given by the Euler angles u. Both the helicity 
s a and the intrinsic parity r\ a of the quark state are specified by the label a — 1,. . . ,4, with 
si,3 = |, S2,4 = — |, ?7i,2 = +1 and 773,4 = —1. The helicity A a of the diquark is specified by 
the label a, with As,o,3 = 0, A + i = +1 and A_i = — 1. The corresponding intrinsic parities 
are 775,3,+!,-! = +1 and r] Q = —1. Note that our use of these labels differs slightly from that 
of Ishii et al. ||18|| . The rotation matrices appearing Eq. ( P7|) are 

' d x l\Q) ^ 
d x l 2 [Q) 



S a a (uj) = exp (is a 'ip) 



exp [ts c 



(28) 



and 



^10 0^ 



= exp(iA a /^) 



1 

d\6) 



exp (i\ 



(29) 



' A„/A a 



where the c?(0) are Wigner ci-functions in Edmonds' convention p7 |. 

These basis states must now be projected onto good angular momentum: 



\p, a, a; JM) 



/2J + 1 



8vr 2 



(30) 



where p = |p|. The resulting Faddeev equation for states of spin J has the form 
XUP0,P) = / dp' p' 2 dp'Ff a b (p ,p,p' ,p')X^ b (p' ,p'). 



(31) 



/9,6 ( 27r ) 

The expression for the kernel is given by the formula (D.9) of Ref. \TE\. We reproduce its 
form in our notation here, 



12 



x 

b.c 



1 

X 



(po + Po) 2 — (p 2 + P' 2 + 2pp' cos 6) — M 2 
x { [ll^i-d)^ + ^rf V2 (-^)7 b 

-m 7 * ^(-^h 6 _ 2d 1 / 2 {-e)d\-e) K , x j 

-2^ 1 / 2 (-^rf 1 (-0) Aa/Ab + 2Mrf 1 / 2 (-^)d 1 (-^) VAb 

+2d l i 2 (-e)d\-e) K , XcPcl b 

+2Yj 1/2 (-9)p' c d\-9) XcXb \S(P/2+p)} , 
xr^P/2-p). (32) 

In Eq. ([32]), the asterisk denotes complex conjugation with respect to the explicit factor of 
i in the spherical-basis components of vectors only. The matrices d l l 2 (6) and d}{6) are the 
ones appearing between the phase factors in Eqs. ( [28]) and ( p9|) respectively. In the spin-| 
channel, the numerical coefficients —3(7*/^ are the ones appearing in Eq. (|25|), with CX 
equal to +1 if both indices refer to scalar diquarks, —1 if both refer to axial diquarks, and 
v^3 if they are mixed. 

Finally the Faddeev equation has to be projected onto positive parity states. The result- 



ing equation has the same form as fl3T|), but with the positive-parity kernel, 



where the quark index /3 is defined such that sp = — sp and m — r]p, and the diquark index 
b is defined similarly. The phase factor is 

zp h = vM-l) J - S ~ jb , (34) 

where T]p and rj^ are the intrinsic parities of the quark and diquark, and their intrinsic spins 
are s = |, j^o = and .73,4.1-1 = 1. The parity projection cuts the set of 20 coupled 
equations down to 10. Of these, two describe only states with spin-| and so decouple from 
the spin- 1 channel to leave 8 coupled equations. 
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A very similar set of equations can be derived for spin-| states. The kernel has the 
same form as given in Eq. (0); only the numerical coefficients differ from spin-| case. From 
Eq. ( p6|) we see that C„/„ is equal to +2 if both indices refer to axial diquarks and zero 
otherwise. Since the components involving scalar diquarks do not contribute, the number 
of coupled equations is again reduced to 8. 



V. NUMERICAL METHOD 

We solve the Faddeev equation (|3T| ) in the rest frame of the nucleon (or A), where 
P = (E,0). To avoid the singularities of the kernel, we perform a Wick rotation on the 
energy variables p and p' : p — * a + ip^ with 

M - m diq 

a = (35) 

where is the lower-energy pole in the diquark T-matrices, (|12j) and (|T6|) (i.e. the mass 
of the lighter diquark) . 

Because of the Wick rotation, we have to solve 8 complex coupled integral equations. 
We do so using the iterative method of Malffiet and Tjon p9| , |30|| . The set of equations may 
be written schematically as 

K{E)<$> = $, (36) 

where the kernel K(E) depends nonlinearly on the energy eigenvalue E. Rather than solve 
this directly, we solve instead the linear eigenvalue problem 

K(E)<£> = \(E)$, (37) 

for a fixed value of E. This is done iteratively, by acting with K{E) on some initial guess for 
the vector $ of vertex functions to obtain a new vector. This is repeated until the vectors of 
functions in successive iterations are simply proportional to each other. The proportionality 
constant is then equal to \{E). We search on E until we find a value for which X(E) = 1 
and so the solution $ satisfies the original equation (p6|) . In the present problem, we use a 
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simple initial guess consisting of a Gaussian for the real part of each of the 8 components 
of $, and the derivative of a Gaussian for each of the imaginary parts. We find that 4 or 5 
iterations are generally enough to determine \{E). 

The parameters of the model that describe the interactions in the qq channels are fixed 
using the solutions to the Faddeev equation at zero density The values of the scalar and 
axial-vector couplings g s and g A that reproduce the masses of the nucleon and the A are 
listed in Table I for our two parameter sets. 

Before describing our calculations at finite density we should compare our zero-density 



results with those in Ref. ||18|| . Note that apart from the difference in cut-off scheme and the 
use of current quark masses, a further difference between our approach and theirs is that 
Ishii et al. do not try to reproduce the nucleon and A masses exactly, but investigate the 
range of parameters g s /gn and g A /g-w that can give bound states with reasonable masses. 
When only the scalar coupling is included (g A / 'g n = 0), we find that the minimum value of 



9s I ' 9ir to get a bound nucleon is 0.8 compared with 0.5 in Ref. UI8"H , both for a quark mass of 
M = 400 MeV. Although these values are rather different, one should note that the nucleon 
is very weakly bound in the absence of the axial qq coupling. For example, Ishii et al. find 
a nucleon binding energy of only about 40 MeV for g s jg-K — 0.8. It is therefore not too 
surprising that the point at which the nucleon becomes bound is rather sensitive to details, 
such as the choice of regulator or the use of the chiral limit. 

With both scalar and axial couplings the nucleon is more strongly bound and one might 
hope that results are less sensitive, but in this case it is harder to make meaningful compar- 



isons with the results of Ref. |18j because of the problem with the erroneous factor of 4 in 
the axial channel mentioned above. For a quark mass of M = 420 MeV we are able to fit the 
nucleon and A masses with g s / 'g n = 0.735 and g A j 'g n = 0.33. The value for the axial- vector 



coupling is smaller than the minimum value to bind the A of 0.44 in Ref. [18|]. The value 
for the scalar coupling is rather somewhat larger than the range considered in that work. 
As in the case of the diquarks, there are indications that our approach tends to give less 
attraction in the scalar channel, compared with that of Ishii et al, but more attraction in 
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the axial channel. Overall, though, the results are qualitatively similar. 

The Faddeev equation at finite density is solved using the same methods. In this case the 
constituent quark mass is density dependent, and the 3- momentum of the valence quarks 
must be restricted to be larger than the Fermi momentum kp. This acts as a lower cut-off 
of the momentum integrals: kp < |pj| < A for the three- momenta pj of all three quarks. A 
similar Wick rotation is performed on the energy variables, but with 

sjM{k F f + k F - E^(k F ) 
a(k F ) = -* , (38) 

where Edi q (kp) is the lower-energy pole in the diquark T-matrices at finite density. 



VI. RESULTS 

We have solved the Faddeev equation at finite density for the energy E N of a nucleon at 
rest. In Figs. [|(a) and |3](b) we show the binding energy of the nucleon as a function of 
the Fermi momentum kp for our two sets of parameters. For these parameter sets the scalar 
diquark is more strongly boundthan the axial, and so the nucleon binding energy is defined 
as 

B N = E s + yfk 2 F + M 2 - E Nl (39) 

with respect to the quark-diquark threshold. This is the relevant threshold since the NJL 
model does not exhibit confinement. We compare Bn to the binding energy of the scalar 
diquark 

B s = 2^Jk 2 F + M 2 — E s . (40) 

The binding with respect to the three-quark threshold is given by the sum of B^ and B$- 
One can see that the behaviors of the nucleon and the scalar diquark at finite density are 
quite different. The binding energy of the diquark is very large and initially tends to increase 
with density, only decreasing when the Fermi momentum approaches the cut-off. In contrast, 
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FIG. 3. The binding energy of the nucleon (continuous curve) and of the scalar diquark (dashed 
curve) as a function of &f for the parameter sets with (a) M = 450 MeV and (b) M = 500 MeV. 
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k F (GeV) 

FIG. 4. The binding energy of the A (continuous curve) and of the axial- vector diquark (dashed 
curved) as a function of kp for the parameter set with M = 450 MeV. 

the binding of the nucleon decreases rather quickly, so that it is only marginally bound at 
nuclear matter density (kp = 270 MeV). 

We have also calculated the binding energy of the A as a function of kp and this is 
shown in Fig. [|. In this case the threshold is determined by the binding energy of the axial 
diquark, which is also shown in Fig. [|. The behavior is similar to that found in the case of 
the nucleon. Since the axial coupling is smaller than the scalar coupling, the axial diquark is 
less strongly bound than the scalar one. Nonetheless it remains bound over the whole range 
of densities considered, while the A becomes unbound near nuclear matter density. Chiral 
symmetry restoration occurs at kp = 359 MeV and kp = 372 MeV for the parameter sets 
with M = 450 MeV and M = 500 MeV respectively, which corresponds to the cusps in the 
diquark binding in Figs. |3| and f|. 

Another perspective can be gained by plotting the total energy of the nucleon instead 
of its binding energy, as is done in Fig. [5]. For comparison, we have also plotted the quark- 
diquark threshold and the three-quark threshold. One can see that the nucleon energy 
increases only slightly with density, while the quark-diquark threshold decreases more and 
more quickly as one approaches the density of chiral restoration (which is basically a con- 
sequence of the vanishing of the constituent quark mass at the transition). The density 
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FIG. 5. The energy of the nucleon (continuous curve) as a function of kp compared with 3/2 
times the mass of the scalar diquark (long-dashed curve) , the quark-diquark threshold (short-dashed 
curve) and the three-quark threshold (dash-dotted curve), for the parameter set with M = 450 
MeV. 

dependence of the A energy shows the same qualitative behavior. 

To answer the question raised in the introduction about the possible competition between 
diquark and three-quark clustering, we have also plotted the quantity \Eg- If En < ^E$ 
then a system of six quarks will prefer to form two nucleons, while if En > ^E$ it will 
form three diquarks. Fig. [5] shows that nucleons are more stable than diquarks only for 
Fermi momenta smaller than about 130 MeV, which corresponds to 1/8 of nuclear matter 
density. For our other parameter set, with M = 500 MeV, nucleons remain more stable than 
diquarks up to kp = 170 MeV, or 1/4 of nuclear matter density. These results clearly cast 
doubt on the validity of the model in this regime. 



VII. CONCLUSIONS 

We find that the NJL model predicts that, except well below nuclear matter density, it is 
energetically much more favorable to form three diquarks than it is to form two nucleons, and 
that there are no nucleon instabilities for the high densities where color superconductivity 
is expected to occur. 
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However we should not conclude on the basis of these results that there is no competition 
between nucleon formation and diquark condensation. Since this model predicts that nuclear 
matter should consist of diquarks, a result clearly at odds with what is actually observed, 
some elements of reality are missing. The most obvious of those is the lack of explicit 
confinement in the NJL model, something that is not easily remedied. It poses a rather 
deep question about the interpretation of the NJL model at finite densities. In much of 
the recent literature it has been used to model an instanton-induced interaction between the 
quarks in high-density matter beyond the chiral/deconfining phase transition. We have tried 
to follow a more traditional approach, where the model is interpreted as a density model 
for hadron structure in the vacuum. In that case the model contains unphysical diquark 
degrees of freedom, which may be ignored at zero density as being "irrelevant". Both of 
these approaches have flaws. In the first (high-density) interpretation, we have no clue as to 
what remnants of confinement may play a role in the interaction. In the second (low-density) 
interpretation we cannot even describe nuclear matter properly. 

Standing back and looking at our results in the light of these problems with the model, 
we note from Fig. ^ that the nucleon energy is roughly independent of density. It is mainly 
the increase in the diquark binding that renders the nucleon unstable. One could naively add 
to the NJL model a three-body force that provides attraction in the color-singlet channel, 
to incorporate an approximate description of confinement. Alternatively one might modify 
our treatment to use confined rather than free quark and diquark propagators [pOR . Both of 
these choices have some appeal, but neither really deals with the underlying mechanism of 
confinement. 

One should also remember that our results have all been obtained using the "rainbow- 
ladder" approximation to the diquark Bethe-Salpeter equation and the Faddeev equation. 
In the context of models with nonlocal interactions between the quarks, it has been shown 



3l| , |32"| that this approximation over-predicts diquark binding. With a toy model for the 
gluon propagator which leads to quark confinement, diquark condensation has been shown 
to occur even though diquarks are no longer bound at zero density |53 . Hence another way 
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to approach this problem may be to study the Faddeev equation within a less restrictive 
calculational scheme. 

In conclusion, it is clear that the question of how confinement and the related three-quark 
correlations affect the phase structure of baryonic matter still remains to be answered. 



ACKNOWLEDGMENTS 



We are very grateful to Noriyoshi Ishii for providing us valuable details of the derivation 
of the Bethe-Salpeter and Faddeev equations. This work was supported by the EPSRC 
under grants GR/J95775 (Manchester) and GR/L22331 (UMIST). 



21 



REFERENCES 

[1] I. Barbour, S. Morrison, E. Klepfish, J. Kogut, M.-P. Lombardo, Nucl. Phys. B, Proc. 
Suppl. 60A, 220 (1998). 



[2] 
[3] 
[4] 
[5] 



[7; 



[9 
[10 

[11 

[12 
[13 
[14 
[15 
[16 

[17; 

[18 
[19 



M. Alford, A. Kapustin and F. Wilczek, Phys. Rev. D 59, 054502 (1999). 

M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B 422, 247 (1998). 

R. Rapp, T. Schafer, E. Shuryak and M. Velkovsky, Phys. Rev. Lett. 81, 53 (1998). 

D. Bailin and A. Love, Phys. Rep. 107, 325 (1984). 

M. Iwasaki and T. Iwado, Phys. Lett. B 350, 163 (1995). 

J. Berges and K. Rajagopal, Nucl. Phys. B538, 215 (1999). 

D. I. Diakonov, H. Forkel and M. Lutz, Phys. Lett. B 373, 147 (1996). 

T. Schwarz, S. Klevansky and G. Papp, Phys. Rev. C 60, 055205 (1999). 

M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. A638, 515c (1998). 

K. Rajagopal, |hep-ph/ 9908350. 



T. Schafer, |nucl-th/99110T7 . 



M. Beyer, W. Schadow, C. Kuhrts and G. Ropke, Phys. Rev. C 60, 034004 (1999). 

A. Buck, R. Alkofer and H. Reinhardt, Phys. Lett. B 286, 29 (1992). 

N. Ishii, W. Bentz and K. Yazaki, Phys. Lett. B 318, 26 (1993). 

H. Meyer, Phys. Lett. B 337, 37 (1994). 

S. Huang and J. Tjon, Phys. Rev. C 49, 1702 (1994). 

N. Ishii, W. Bentz and K. Yazaki, Nucl. Phys. A587, 617 (1995). 

C. Hanhart and S. Krewald, Phys. Lett. B 344, 55 (1995). 

22 



[20] G. Hellstern, R. Alkofer, M. Oettel and H. Reinhardt, Nucl. Phys. A627, 679 (1997); 
M. Oettel, G. Hellstern, R. Alkofer and H. Reinhardt, Phys. Rev. C 58, 2459 (1998). 



[21 
[22 
[23 
[24 
[25 
[26 
[27 

[28 
[29 

[30; 

[31 
[32 
[33 



M. C. Birse, J. A. McGovern, S. Pepin and N. R. Walet, |nucl-th/99U5U3l 



Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961); 124, 246 (1961). 

S. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 

T. Hatsuda and T. Kunihiro, Phys. Rep. 247, 221 (1994). 

C. Itzykson and J.-B. Zuber, Quantum Field Theory, (McGraw-Hill, New York, 1980). 
N. Ishii, private communication. 

A. R. Edmonds, Angular Momentum in Quantum Mechanics, (Oxford University Press, 
Oxford, 1957). 

M. Jacob and G. C. Wick, Ann. Phys. 7, 404 (1959). 

R. A. Malfliet and J. A. Tjon, Nucl. Phys. A127, 161 (1969). 

G. Rupp and J. A. Tjon, Phys. Rev. C 37, 1729 (1988). 

A. Bender, C. D. Roberts and L. van Smekal, Phys. Lett. B 380, 7 (1996). 

G. Hellstern, R. Alkofer and H. Reinhardt, Nucl. Phys. A625, 697 (1997). 

J. C. R. Bloch, C. D. Roberts and S. M. Schmidt, |nucl-th/9907086 . 



23 



